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Abstract 

The adsorption at the interface between an aqueous solution of several surface- 
active agents and another fluid (air or oil) phase is addressed theoretically. We 
derive the kinetic equations from a variation of the interfacial free energy, solve 
them numerically and provide an analytic solution for the simple case of a linear 
adsorption isotherm. Calculating asymptotic solutions analytically, we find the 
characteristic time scales of the adsorption process and observe the behavior of the 
system at various temporal stages. In particular, we relate the kinetic behavior of 
the mixture to the properties of its individual constituents and find good agreement 
with experiments. In the case of kinetically limited adsorption, the mixture kinetics 
is found to be considerably different from that of the single-surfactant solutions 
because of strong coupling between the species. 



To be published in Langmuir. 
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1 Introduction 



The kinetics of surfactant adsorption plays an important role in various interfacial phe- 
nomena and has been an active field of research, both experimentally and theoretically, 
since the 40s [0. Recently, a new theoretical approach based on a free-energy formalism 
was presented and applied to nonionic as well as ionic surfactant adsorption |J . The main 
advantage of this approach is that all the equations are derived from a single functional. 
This feature facilitates generalizations of the model to more complicated situations. 

Surfactant mixtures are used in numerous industrial applications, and are also encoun- 
tered in many systems because of the presence of surface-active impurities. Experiments 
portray a large variety of phenomena specific to the kinetics of mixed systems For 
instance, more complex dynamic surface tension is observed due to competition between 
the different species. 

The equilibrium behavior of mixed-surfactant solutions was studied in detail in pre- 
vious works P|-||13||. One of the important results, both theoretically and from the ap- 
plication point of view, is the ability to relate the mixed-surfactant behavior to that of 
the better understood, single-surfactant one. In the current work we focus on kinetic 
aspects, deriving equilibrium results merely for completeness and comparison with previ- 
ous models. One of our aims is to predict the mixture kinetics from the behavior of the 
single surfactants. A particularly interesting question is whether mixing several species 
would lead in certain cases to a significant difference in the kinetics as compared to the 
single-surfactant systems. 

In this paper we extend our previously introduced model to describe the competition 
between two nonionic adsorbing species. In sec. |2| we present the model and derive the 
equilibrium relations and kinetic equations. The complete set of kinetic equations can be 
solved only numerically, as is done in sec. |3|. We then discuss, in sec. f|, the asymptotic time 
dependence of diffusion-limited and kinetically limited adsorption and the corresponding 
characteristic time scales. We focus on the relation between the adsorption behavior of 
the mixture and the properties of its individual constituents. Good agreement is found 
between the experimental results and our predictions. In addition, we give an analytic 
solution for the kinetics in the simple case of a linear adsorption isotherm. 

In order to make the formulation as transparent as possible, we have tried to minimize 
the number of symbols and used dimensionless forms whenever possible. To facilitate 
translation into experimentally useful, dimensional quantities, a nomenclature table is 
provided in Appendix B. 

2 The Model 

Our model is based on the free-energy formalism presented in detail in previous papers 
0. We consider a semi- infinite aqueous solution of nonionic surfactants having, at x — 0, 
a flat, sharp interface with a non-polar fluid phase (such as air or oil). The solution is in 
contact with a bulk reservoir, at x — > oo, containing two types of surfactant molecules, de- 
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noted by A and B. It is considered dilute, with the volume fraction of the two constituents 



well below their critical micelle concentration (cmc) values [14]. At the interface, how- 



ever, the surfactant volume fraction is usually much larger. Hereafter we use <fi to denote 
volume fraction of surfactant A, ip for volume fraction of surfactant B, /i for chemical 
potential and D for diffusion constant. The subscripts A and B distinguish between the 
two different surfactants. The subscripts 0, l,b are used to distinguish between different 
positions in the solution, corresponding, respectively, to the interface, sub-interface layer 
and bulk reservoir. 

The excess interfacial free energy is equal to the reduction in surface tension, 

A 7 [0,^]=/ {A/[0(x)] + A/[^(x)] }dx + / o (0 o ,^o) (2.1) 



where Af is the bulk contribution per unit volume of each species, and f the interfacial 
contribution per unit area. In dilute nonionic surfactant systems, the dominant contribu- 
tion to the surface tension is usually the interfacial one, fa. Since the surfactant solution 
is considered dilute, the bulk free energy is taken to be the sum of A _/((/>) and Af(i(>). 
Each of these contributions contains the entropy of mixing in the ideal, dilute limit and 
a contact with the bulk reservoir. For species A we write 

A/(0) = — [4> In cf> - 4> - (0 b In b - 4> h ) - // bjA (0 - b )] (2.2) 

and similarly for species B. As this is a lattice model, it is conveniently formulated using 
volume fractions as the degrees of freedom. The translation to mole fractions or concen- 
trations requires specifying molecular sizes for the various species. Here, the surfactant 
molecular size, a, is assumed to have the same value for both species, on account of sim- 
plicity ||31)|| . All the energies, chemical potentials and interaction parameters are expressed 
in units of k^T where &b is the Boltzmann constant and T the temperature. At the in- 
terface, the surfactant volume fractions are usually much larger and we have to consider 
the full, non-ideal entropy of mixing and additional interaction terms, 

/o(0o, ^o) = -rr [0o In 0o + V'o hi ipo + rj Q In r] - (a A + ^i,a)0o- ("b + A*i,b)^o 



Pa , 2 Pb i 2 ii 



(2.3) 



where a is the energy gain of adsorption onto the interface, P the interaction energy 
between molecules of the same kind, e the interaction between different surfactants fT6|| , 



and TjQ = 1 — 4>q — ipQ is the surface coverage of the solvent (water). The interface is in 
contact with the sub-surface layer, having a chemical potential /^i = fi(x — > 0). Out of 
equilibrium fii may differ from the bulk chemical potential, /ib- Although both a and fii 
are linearly coupled with the surface coverage, their physical origin is quite different. The 
former is constant in time, characterizing the surface activity of the specific surfactant 
and mainly dependent on the molecular structure (number of hydrocarbon groups in the 
surfactant tail). The latter is a time dependent function participating in the surface 
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kinetics. Variation of A7 with respect to and ip yields the excess chemical potential at 
distance x from the interface, 



A/i A (x) = ln<f)(x) -/i b ,A 

A/j, b (x) = lnV(x) -Mb,B (2.4) 



and at the interface itself, 



A/i ,A = In — - a A - /?A0O - #0 - A*1,A 
^0 

A/io.B = In — - a B - A^o - - A*i,b (2-5) 
Vo 

2.1 Equilibrium 

At equilibrium the chemical potentials are equal to their bulk values throughout the 
solution, leading to two uniform profiles, (j)(x > 0) = 0b and ip(x > 0) = iph- At the 
interface we obtain a Frumkin adsorption isotherm [|17]], generalized for the A/B mixture 

case: 



, = 0b(l ~ ^0) 

W ° ^ b _|_ e -(a B +/3 B V'0+^o) ^ ' 

The adsorption of species A depends on the species B because of the entropy of mix- 
ing (steric effect) and surfactant-surfactant interactions. The corresponding generalized 
Langmuir isotherm is obtained in the limit of no interactions, /5a = Pb — £ — 0. Finally, 
the equilibrium equation of state, A7 = A 7 (0, -0), takes the form 

A 7 = i (lnr, + ^0 O 2 + y^o 2 + e<M^ (2.7) 

which is equivalent to Lucassen-Reynders' result ||, when differences in molecular sizes 
are neglected. 

2.2 Out of Equilibrium 

We apply the procedure presented in Ref . [Q] , where the kinetic equations are derived from 
the variation of the free energy. The procedure is generalized for the present case of a two 
nonionic surfactant mixture. Since the bulk solution is dilute, two independent diffusion 
equations for the two surfactants are obtained, leading to two Ward-Tordai equations |T3|] , 
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similar to previous models |Q: 




+ 



+ 2^ b - Vi 



(2i 



At the interface, however, the two species are correlated and the procedure yields two 
coupled kinetic equations: 



<90o 
dt 


= 01—5- 

a z 


In 


^01 

V0o 


dijj 
dt 


= V'i— 

er 


In 


^i 



-rVo] + "a + /5a0o + #o 



(2.9) 



As can be seen from the equation above, the coupling between the kinetics of the two 
species arises from an interaction term as well as an entropic one (via the 770 term). It 
should be mentioned that this form of interfacial kinetic equations is particular to our 
free-energy approach. Close enough to equilibrium, however, it coincides with the usual 
adsorption-desorption form, as discussed in ref. |§. The system of four equations, (|2.8|) 
and (|2.9|) , with the appropriate initial conditions, completely determines the mixture 
kinetics and equilibrium state. Several limits, such as diffusion-limited and kinetically 
limited adsorption, can be treated analytically, as presented in sec. 4. The full solution 
of the mixed-kinetics case is obtained numerically. 



3 Numerical Solution of the Full Equations 

Several numerical schemes have been proposed for solving the Ward-Tordai equation with 
various boundary conditions [fHfl-fZT]. We generalized the recursive scheme suggested by 



Miller et al. |19j to a surfactant mixture having time dependent boundary conditions. The 
complete set of integro-differential equations were solved, as is explained in the Appendix. 

The time dependence of 0o, tpo and their sum 0o + ipo can be seen in Fig. 1. The 
mixture parameters are specifically chosen to show the interesting case of competition 
between the two species. While surfactant B diffuses more rapidly and is more abundant 
at the interface during the initial stages of the adsorption process, surfactant A has 
a higher surface affinity and dominates the later stages of adsorption. We note that, 
because of this competition, surfactant A not only takes over the adsorption at the later 
time stages, but it also forces surfactant B to desorb from the interface. In Fig. 2 the 
dynamic surface tension is shown for the same time scales and adsorption parameters as 
of Fig. 1. The competition between the A and B surfactants results in a more complex 
decrease of the surface tension at intermediate times. 
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4 Limiting Time Behavior 



The asymptotic time behavior close to equilibrium of the kinetic equations provides a 
means for comparing theoretical predictions to experimental results ||. Two limiting 
cases of adsorption can be distinguished and treated separately for long and short times: 
Diffusion- Limited Adsorption (DLA) occurs when the kinetics at the interface is much 
faster than the diffusion from the bulk to the sub-surface layer. Equation ( |2.9| ) then 
equals zero 0, since the interface is taken to be at equilibrium with the sub-surface layer. 
The profiles in the bulk, 0(x > 0) and ip{x > 0), depend on the diffusion and evolve 
in time. On the other hand, Kinetically Limited Adsorption (KLA) takes place in the 
opposite limit, for which <f>i = (j) h and ipi = ip^ throughout the process, implying that the 
bulk is constantly at equilibrium with the reservoir but not with the interface. 



4.1 DLA process at long times 

Diffusion-limited processes such as the DLA considered here have a characteristic asymp- 
totic £ -1 / 2 dependence as t goes to infinity. This is demonstrated in Fig. 3. The asymptotic 
time dependence of all variables can be written in a generic form as 



0q(*) - 0o, eq (l - \froft\ 



(t) ~ b [l-Jn/t 



A 7 (f) ~ A 7eq (1 - y/r^/t) (4.1) 

where the subscript 'eq' stands for equilibrium values. Note that the various time con- 
stants may generally depend on the surface coverages of both species. Taking the limit 
t — > oo of the Ward-Tordai relation (|2.8|) we get 




ri,A~ ^ (4.2) 



The presence of surfactant B changes the equilibrium value of O [ see ec i- Q2.6Q 1. This is 
the only reason why the relaxation time t^a depends on surfactant B. For a DLA process, 
the surface and sub-surface layers are considered at equilibrium with each other, and the 
coverage instantaneously follows any change in the sub-surface concentration according 
to 

0i(l-V>o) 



^l(l - 0p) , V 

^° fa + e -(a B +Mo+s<t> ) \ 6 > 
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We now omit the subscript 'eq' for brevity. Substituting ( |4.1| ) into ( |4.3| ) and comparing 
powers of t, we get two linear equations for t 0) a and r B 



Voy/nJ = (1 - ipo - Pa4>oVo) v^a + i>o (1 - ^7o^) v 7 ^" 

^Ov 7 ^ = (1 - 00 - Pb^oVo) y/T^B + 00 (1 ~ ^) y/TQ~A (4.4) 

Let us consider a single- surfactant system (say the A) as a reference state. This is simply 
achieved by setting ^o.eq and t 0i b to zero, and denoting O - ► 0o, Vo — ► Vo = 1 — 0o- 
Consequently, 

1 - ~ 



f ° A = i — — ITTV fl ' A ^ 4 - 5 ) 

1 - /3 A (1 - 0oJ 0o 

where the "bars" throughout the paper denote values obtained for this single-component 
system. 

In a DLA process, the equation of state ( |2.7|) holds also out of equilibrium p|, since 
the equilibration of the interface with the sub-surface layer is very fast for both species. 



Substituting ( |4.1|) into (|2.7| ), the asymptotic behavior of the surface tension is obtained: 



- « 2 A7\A^ = ^ - /3 A 0o - £0oV>oJ y^A + ~ Mo - £<t>oi>oJ s/t^b (4.6) 

Without surface interactions (Pa = Pb = £ = 0) only steric effects and surface ac- 
tivity terms are taken into account. The Frumkin-like isotherm (|4.3|) is simplified in this 
approximation to a Langmuir-like isotherm, and eq. (f4.4j) takes the form 



\/^a = (1 - 0o ) \/n^A - ^o^/t^b 

v /r^i = (1 - tp ) V /T^i - Ov /tXa (4.7) 

Using eqs. ( |4.5|) -( |4~?7|) we obtain a simple expression, relating r 7 of the mixture with those 
of each species separately, t 7i a and f 7jB , 



^7x7^7 = A 7a(0o/0o) \/r 7 ,A + A7b(V>oM)) \/ t 7,b (4. 



where A7a and A7b are the equilibrium reduction in interfacial tension of the single- 
surfactant solutions. This correspondence relates the time scale of the surface tension 
relaxation in the mixture with the time scales corresponding to the individual species. 
As was pointed out previously, most common nonionic surfactants usually undergo DLA. 
Thus, the above result provides a convenient way of predicting the behavior of multi- 
component surfactant mixtures based on single- surfactant data. In Table 1, we compare 
the predicted y/f^ of eq. ( |4.8|) with experimental results obtained by Fainerman and 
Miller [|J for a sequence of Triton X mixtures. Based on single-surfactant values and 
equilibrium isotherms for the mixture, the two terms of eq. ( [4.8| ) are calculated separately. 
The agreement between theory and experiment is quite good, although experiments were 
limited to cases having one species dominating the adsorption. The last entry in the 
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table corresponds to a mixture of Triton X-405 and Triton X-165. Here the predicted r 7 
deviates from the experimental one by 33%. Equilibrium measurements on this mixture 
reveal an increase in the X-165 coverage in the presence of the X-405 |4j], implying strong 
interfacial interactions between the species. The deviation in the predicted kinetics in 
Table 1 probably arise from those interactions, which are not taken into account by 
eq. ( |4.8|) . It is possible to treat also the general case, including interactions between 
surfactants, by using the full equations (|4.4j )-( [L6|) instead of the simplified one ( |4.8|) . 
Such a procedure, however, involves three additional fitting parameters (/?a, Pb and e). 
Nevertheless, as demonstraed in Table 1, the simple prediction ( [4.8|) may be applicable 
to various experimental systems. 



4.2 KLA process at long times 

Although most nonionic surfactants undergo a DLA process, the adsorption of some 
nonionic surfactants as well as ionic ones is found to be kinetically limited. For a KLA 
process, the bulk of the mixture is assumed to be at equilibrium, 0(x) = 0b and ip(x) = -0b- 
The equations governing the kinetics are now the interfacial ones, (|2.9|) . 

Asymptotic solutions at t — ► oo of first-order equations such as ( |2.9|) have an exponen- 
tial form characteristic of KLA. Linearizing eqs. ( |2.9|) about the equilibrium state, 0o, e q 
and f/Veq, two time scales denoted r+ and r_ emerge (r_ > t+). These collective time 
scales correspond to the kinetics of a certain combination of surfactant coverages, 

CiA^o + CaA^o ~ e- t/T - 

C 3 A0o + C 4 A^o ~ e^ /T + (4.9) 

where A0 O = O — 0o, eq , A-0 O = ip — V>o,eq, and C\ . . . C4 are constants. Since r_ > r + , it 
is t_ which limits the kinetics of the system. 

In the simple case of no surface interactions between A and B ( /3a = Pb = £ = 0, yet 
keeping steric and surface activity effects), the expressions for r± are 

2_ _ 1 - ^0 | l-0o 
where ta and tb are the KLA time 

tb 

yet 0o and tpo are the equilibrium values for the mixture. The behavior of the mixed 
system combines the single-surfactant kinetics in a complicated manner. We can gain 
some insight on this coupling by considering some simple cases. 



| -^o + Wo\ _ 4 ^p_ (410) 



ta t b I t a t b 



scales of the single surfactant case, defined as 



a 2 ( 0o \ 2 



Da \4>b, 

a 2 (jjQ 
Db Ub, 



g-«A 



e-° B (4.11) 
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For the case where the interfacial kinetics of surfactant A is much slower than that of 
B, ta 3> Tb, eqs. (|4.9| ) and ([4.10|) are simplified to 



(1 - 0O,eq)A0 O - V'CeqAV'O ~ e * /r ~ \ T_ = T A (l - / T]q 

A^ ~e- t/T+ ; r + = r B /(l-0 o ) (4.12) 
In the other limiting case, where the two species have similar time scales, ta — tb, we get 

A0 O - A^o ~ e~ t/r - ; r_ = r A /r] 

0O,eqA0 O + V'O.eqA^o ~ e"*^ 5 T+ = T A (4.13) 

The factor l/r/o in r_ is quite interesting. Since the equilibrium surface coverage of 
the solvent, 770, is usually very small in surfactant systems, this factor implies that the 
coupling in a surfactant mixture undergoing KLA may lead to a significant reduction in 
adsorption rate. In this regime the mixture behavior may differ considerably from that of 
its individual constituents. Because of the relatively large factor of 1/t]o, the KLA time 
scale may exceed the DLA one and the adsorption would become kinetically limited. 



4.3 Short time behavior 

DLA regime. Here we assume that the surface and sub-surface layers are already at 
equilibrium and examine the system at time scales short compared with the diffusion 
mechanism. Note that the DLA behavior cannot strictly start at £ = 0, since at that 
instance the interface and sub-surface layers are not at equilibrium with each other. As- 



suming a DLA time dependence of the form 4>o(t) ~ const. + yt/r A , the const, is found 
to be roughly equal to 20b. In other words, only once the surface coverage has exceeded 
a value of 20b, indicating an almost complete depletion of the sub-surface layer fi~8| , 



can one assume a process limited by diffusion (the same argument applies to the second 
surfactant). Figure 4 shows the surface coverage and surface tension at this DLA stage 
plotted as function of t 1 ^ 2 . The linear dependence of the surface coverages is evident, 
whereas the dynamic surface tension, having contributions from both species, exhibits a 
small deviation from the t 1//2 behavior. Since surfactant A has been taken to be much 
slower than surfactant B, the t 1//2 behavior of the two surfactants overlaps only for a very 
short period of time. 

Beyond the DLA regime. At the very early time stages of the process, most of the 
molecules in the sub-surface layer migrate rapidly to the interface. Only when the sub- 
surface layer becomes nearly depleted, do molecules from the bulk start migrating towards 
the interface by a diffusive mechanism. Before diffusion sets in, the DLA assumption is 
invalid and the interfacial kinetics must be considered explicitly. To address these very 
early time stages we assume that the bulk solution is still at its initial equilibrium state, 
unperturbed by the presence of the interface. The leading time behavior of the surface 



9 



coverage is found from ( |2.9| ) to be linear, as can be seen in Fig. 5, 

<f) (t -> 0) ~ b + —IrfoaAt 
a 2 

(t -> 0) ~ ^ b + ^ Vb« B * (4. 14) 

or 

Any correlation between the two species vanishes during this initial stage since the sur- 
factants at the interface are dilute and non-interacting. We would like to emphasize that 
those very early time scales (10 _11 -10 -7 s) are of no experimental interest. We study 
them as an analytic limit of the general kinetic equations. 



4.4 The linear isotherm limit: weakly adsorbing surfactants 

We present now a very simplified and restricted case which, nevertheless, is of interest 
since it has an analytic solution. The assumptions for which this case can be treated are: 
(i) one of the surfactants (say A) has a low surface activity and undergoes a slow DLA 
process; (ii) the second surfactant (B) is assumed to be at equilibrium; and (iii) the two 
species do not interact. Since the surface coverage of surfactant A is always very low, its 
isotherm can be taken as linear in 0! (Henry's law): 

0o = e QA 0i (4.15) 

as can be seen from taking the low-coverage limit of eq. ( |2.6| ). With the above assumptions, 
the time-dependent surface coverages of the two species can be explicitly calculated 



(f>o(t) = [1 -tanh(a A /2)e' /r erfc( v /t7^)]0 b e aA 

Mt) = M-MH (4.16) 

ipb + e 

where r = (a 2 /D^)(l + e aA ) 2 . Note that in the limit t — > 0, the DLA process starts at 
<po(t = 0) = 20b/ (e _QA + 1) which is somewhat smaller than 20b. This result is consistent 
with an early-stage dynamics, where the sub-surface layer is entirely depleted, as was 
discussed in Sec. 4.3. Furthermore, the linear relation ( |4.15| ) is always valid for short 
enough times, when the surface coverage is low (with no further limiting assumptions), 
as demonstrated in Fig. 6. 

5 Conclusions 

We have examined the behavior of nonionic surfactant mixtures by solving numerically the 
governing kinetic equations, as well as deriving analytical limiting cases. A comprehen- 
sive description of the different adsorption stages is obtained. For mixtures of common, 
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nonionic surfactants, the adsorption process can be divided into four temporal stages. 
At the very early times of the adsorption, the process begins with a short stage, where 
the surface coverage and surface tension change linearly with time because of interfacial 
kinetics (Fig. 5). As discussed is sec. 4.3, this stage is practically too short to be observed 
experimentally (usually less than microseconds). During this early stage, the sub-surface 
layer becomes nearly empty, which in turn drives a second, diffusion-limited stage, where 
molecules of both species diffuse from the bulk with a t 1//2 time dependence (Fig. 4). Dur- 
ing this second stage (discussed above in sec. 4.3), the coverage is dominated by the more 
mobile species having a faster diffusion. In cases where the less mobile species is more 
surface active, a third stage is predicted. Here one species undergoes desorption, while the 
coverage gradually becomes dominated by the second, energetically favorable surfactant. 
This competition can be seen in Fig. 1. The final relaxation towards equilibrium is usually 
diffusion-limited, exhibiting an asymptotic t _1//2 behavior (Fig. 3), which was explained 
in sec. 4.1. 

For surfactant mixtures exhibiting kinetically limited adsorption, we find a "syner- 
gistic" effect, where the mixture kinetics may be considerably different from that of the 
individual species. In cases of high equilibrium surface coverage, our model predicts a sig- 
nificant increase in the limiting time scale due to coupling between the two surfactants. 

We have managed to relate the kinetic behavior of the two-surfactant mixture with the 
properties of its individual constituents, [eq. fl4.8|) 1. The kinetic behavior of the mixture 
can be predicted based on equilibrium isotherms and single- surfactant data. Our results 
are in good agreement with experiments. However, since previous experiments were re- 
stricted to surfactants of very different adsorption time scales, our theory could be checked 
only in those experimental conditions. Further experiments on surfactant mixtures, es- 
pecially with comparable adsorption time scales and equilibrium surface coverages, are 
needed in order to get a better verification of our model. 

The problem we have dealt with in this work is another example of the ease at which 
the free-energy approach to adsorption kinetics can be generalized. Examples for further 
extensions are, for instance, the addition of adsorption barriers and treatment of surfactant 
adsorption from micellar solutions. 
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Nomenclature 



Below is a summary of the symbols used in this work. In order to make the formulation as 
concise as possible, we have extensively used dimensionless forms. Whenever necessary, a 
translation to the more practical, dimensional quantities is provided, using the symbols 
of e.g. ref. |§. 

Symbols 

a molecular dimension (equal to the inverse square root of the maximum 

surface density, T^ 1 / 2 ) 
D — diffusion coefficient 

Af — excess free energy per unit volume of the bulk solution 
/o — interfacial contribution to the free energy per unit area 
t — time 

x — distance from the interface 

a energy (in units of k^T) gained by a surfactant molecule by migrating to 

the interface 

(3 energy (in units of k^T) of lateral attraction between two surfactant 

molecules of the same species 
A7 change in interfacial tension 

e energy (in units of k&T) of lateral attraction between two surfactant 

molecules of different species 
\i — chemical potential (in units of k&T) 
t — time scale of a kinetic process 

<ft, if) volume fractions of the two surfactant species (equal to the corresponding 

concentrations multiplied by a 3 ; = a 3 c A , if) = a 3 ce) 

4>o,ipo surface coverages of the two surfactant species (equal to the corresponding 

surface densities multiplied by a 2 ; O — cl 2 T^, ipo = a 2 r B ) 

00,^0 — surface coverages in the single-surfactant solutions 

Subscripts 

A,B — corresponding to surfactant A, B 

— value at the interface 

1 — value at the sub-surface layer of the solution 
b — value at the bulk reservoir 

eq — value at equilibrium 
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Appendix: Numerical Scheme 



The objective of the numerical scheme is to solve the set of four equations — the two 
Ward-Tordai equations Q2.8|) and the two interfacial ones ( [2.9|) . The hardcore of the 
problem is the evaluation of the integral in eq. Q2.8|) , 



1 D A 



ay it 



20 b Vt - f 
Jo 



Mr) 



dr 



+ 



Vt^T 

This is done by discretizing the time variable, t, and transforming the integral into a 
recursion relation JT9|. The singularity in the integrand is removed by substituting x = 
y/t — r. Evaluating the integral at discrete time steps, At = t k+ i — t k , we obtain 



>1 T 



-At -> [01 (*n - tk) + h(tn ~ tk-i)] (Vh ~ Jt k - 



k=l 



where t k = kAt and t n = t. The initial conditions are <f>o(0) = 0i(O) = 0& and V'o(O) — 
ipi(0) = ipb- Adding the rest of the terms in the Ward-Tordai equations, we obtain linear 
relations between O an d 0i, and between ip and ipi. Substituting these relations into the 
kinetic boundary conditions ( p.9|) and using a finite difference to evaluate the first time 
derivative, dcpo/dt — > [0o(tfc+i) ~ 0o(^fc)]/At, a set of four algebraic equations is obtained, 
relating the values of O , if) , 0i and -01 at the current time step, t = t n , 



00 (*n) 



1 D t 



a V vr 



n— 1 



+ : 



20b V^n ~~ X! 01 ~~ tk){ytk+l ~ \ftk- 
1 



k=l 



<n-l 



a V vr 



00 (tn) 



(t n _i) + ^At<j0!(t n )ln 



00 (tn) 



+ a A + /?A0o(^n) + £^o(^n) 



and similar two equations for surfactant B. This set of equations is solved numerically, 
using a binary search method. 
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Table Caption 

Table 1 Comparison of the predicted ^Jr^ [eq. (|4.8|) 1 with experimental results [|J. The 
materials used were sequences of Triton X mixtures. The single- surfactant values, 
0o, ipo, A^av^a") A7bV^b, and equilibrium isotherms for the mixture, 0o, eq and 
-0o,eq, were taken from the same reference. The values for A'-f^/r^ (given in units of 
dyne sec^/cm) are obtained experimentally from the asymptotic slope of the 7 vs. 
t^ 1 / 2 curves [see eq. The predicted values for A7 v /t^ of the mixture and the 

corresponding experimental results are given in the columns indicated by 'th' and 
'exp', respectively. The last column shows the respective error between theory and 
experiment. 

Figure Captions 

Fig. 1 Surface coverage in a mixture of interacting surfactants. The dotted, dashed and 
solid lines are the surface coverages of surfactants A (</>o), B (-0 O ) and the total 
coverage (<f> + ipo), respectively. The assigned parameters are: 0b = 1 x 10~ 4 , 
4 = 2x 10" 4 , a A = 10, a B = 9, f3 A = Pb = 3, e = 1, y/Dl/a = 300 sec~ 1/2 , 
V-Db/a = 900 sec _1,/2 . This implies that surfactant A diffuses more slowly but is 
more surface active (a a > «b)- 

Fig. 2 Dynamic surface tension for the system of fig. 1. 

Fig. 3 Dynamic surface tension for the system of fig. 1, redrawn for long times. The curve 
exhibits the t~ 1 / 2 asymptotic behavior characteristic of DLA. 

Fig. 4 Surface coverage (a) and dynamic surface tension (b) for the system of fig. 1, re- 
drawn for intermediate short times. The curve exhibits the expected t 1 / 2 behavior 
characteristic of DLA at short times. 

Fig. 5 Dynamic surface tension for the system of fig. 1, redrawn for extremely short times. 
The curve exhibits the predicted linear behavior. 

Fig. 6 Comparison between the numerical solution of the complete set of equations (solid 
line) and the analytic solution assuming a linear adsorption isotherm (dashed line). 
The parameters used in the calculations are: 0b = 2 x 1CT 5 , -0b = 10~ 4 , «a = 6, 
a B = 5, p A = p B = e = 0, yTDija = 2 x 10 4 sec" 1 / 2 and y/D^/a = 10 4 sec" 1 / 2 . The 
linear approximation is valid only in the beginning of the process, and fails when 
the coverage becomes large. 
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